attorney Docket No.: PD-201200 
Customer No.: 020991 



Patent 



TITLE OF THE INVENTION 

METHOD, SYSTEM, DEVICE AND COMPUTER PROGRAM PRODUCT FOR A 
DEMODULATOR WITH COMMUNICATIONS LINK ADAPTATION 

CROSS REFERENCE TO RELATED DOCUMENTS 
[01] The present invention is related to commonly owned United States Patent Application 
Serial Number 09/978,291 of Liang et al, entitled "METHOD, DEVICE AND COMPUTER 
PROGRAM PRODUCT FOR A DEMODULATOR USING A FUZZY ADAPTIVE FILTER 
(FAF) AND DECISION FEEDBACK," filed on October 16, 2001 and includes use of various 
technologies described in the references identified in the appended LIST OF REFERENCES 
and cross-referenced throughout the specification by numerals in brackets corresponding to 
the respective references, the entire contents of all of which are incorporated by reference 
herein. 

BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION 

[02] The present invention generally relates to satellite communications systems and more 
particularly to a method, system, device and computer program product for a demodulator 
with communications link adaptation. The present invention includes use of various 
technologies described in the references identified in the appended LIST OF REFERENCES 
and cross-referenced throughout the specification by numerals in brackets corresponding to 
the respective references, the entire contents of all of which are incorporated herein by 
reference. 

DISCUSSION OF THE BACKGROUND 

[03] In recent years, with the development of the third-generation (3G) wireless networks, 
General Packet Radio Services (GPRS) [7] wireless networks are being implemented to 
support multimedia personal communication services. With these services, people are able to 
use personal computers (PCs), laptop PCs, personal digital assistants (PDAs), personal 
information assistants (PIAs), etc., to access multimedia information anywhere. GPRS 
wireless networks use a short burst format to reduce end-to-end transmission delay. Kim and 
Cox [4] proposed a dual mode blind equalizer based on a Constant Modulus Algorithm and 
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their scheme is applicable to short burst transmission formats. Viterbi and Viterbi [6] 
proposed a semi-blind demodulation algorithm, but training sequences are needed to remove 
the ambiguity of inverse tangent function and the algorithm is typically applicable for only fri- 
ary Phase-Shift Keying (M-PSK). 

f 04] However, the above techniques typically are not applicable to a plurality of 
modulation techniques, such as M-PSK, Quadrature Phase Shift Keying (QPSK), Quadrature 
Amplitude Modulation (QAM), Pulse Amplitude Modulation (PAM), etc., and typically treat 
link adaptation and demodulation as two separate tasks. 

[05J Therefore, there is a need for a method, system, device and computer program product 
for a demodulator applicable to M-PSK, QPSK, QAM, PAM, etc., and that combines link 
adaptation and demodulation into one framework. 

SUMMARY OF THE INVENTION 

[06] The above and other needs are addressed by the present invention, which provides an 
improved device, system and computer program product for a demodulator with 
communications link adaptation that is applicable to a plurality of modulation techniques, - 
such as M-ary Phase-Shift Keying (M-PSK), Quadrature Phase Shift Keying (QPSK), 
Quadrature Amplitude Modulation (QAM), Pulse Amplitude Modulation (PAM), etc., as 
compared to conventional demodulators. 

[07] Accordingly, in one aspect of the present invention there is provided an improved 
method, system, device and computer program product for a demodulator with- 
communications link adaptation, including receiving a modulated signal over the 
communications channel; extracting clusters from the modulated signal based on an 
unsupervised clustering technique; computing a mean and standard deviation for each 
extracted cluster; determining categories for each extracted cluster based on a training 
sequence included in the modulated signal; and demodulating the modulated signal based on 
the mean, the standard deviation and the determined categories. 

[08J Still other aspects, features, and advantages of the present invention are readily 
apparent from the following detailed description, simply by illustrating a number of particular 
embodiments and implementations, including the best mode contemplated for carrying out 
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the present invention. The present invention is also capable of other and different 
embodiments, and its several details can be modified in various respects, all without 
departing from the spirit and scope of the present invention. Accordingly, the drawing and 
description are to be regarded as illustrative in nature, and not as restrictive. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[ 09] The present invention is illustrated by way of example, and not by way of limitation, 
in the figures of the accompanying drawings and in which like reference numerals refer to 
similar elements and in which: 

[1 0] Figure 1 is a system diagram illustrating an exemplary satellite communications 
system, which may employ a demodulator with communications link adaptation, according to 
the present invention; 

[1 1] Figure 2 is a block diagram illustrating a demodulator with communications link 
adaptation, which may be used in the system of Figure 1, according to the present invention; 

[1 2] Figure 3 is a block diagram of a system model used to evaluate the performance of the 
demodulator with communications link adaptation of Figure 2, according to the present 
invention; 

[13] Figure 4 is a diagram illustrating a burst format used in GPRS wireless networks, 
according to the present invention; 

[14] Figure 5 is a flow chart illustrating the operation of the demodulator with 
communications link adaptation of Figure 2, according to the present invention; 

[15] Figures 6a and 6b are graphs illustrating the performance of the demodulator with 
communications link adaptation of Figure 2, (a) RMSE in terms of mean and (b) RMSE in 
terms of std, according to the present invention; 

[16] Figure 7 is a graph illustrating the performance of the demodulator with 
communications link adaptation of Figure 2 with respect to error of SQI estimation, according 
to the present invention; 
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[17] Figure 8 is a graph illustrating the performance of the demodulator with 
communications link adaptation of Figure 2 and a BPE-based demodulator when the Rician 
factor K = 9dB widfd=lOHz, according to the present invention; and 

[ 1 8] Figure 9 is an exemplary computer system, which may be programmed to perform one 
or more of the processes of the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[19] A device, method and computer program product for a demodulator with 
communications link adaptation, are described. In the following description, for purposes of 
explanation, numerous specific details are set forth in order to provide a thorough 
understanding of the present invention. It is apparent to one skilled in the art, however, that 
the present invention may be practiced without these specific details or with an equivalent 
arrangement. In some instances, well-known structures and devices are shown in block 
diagram form in order to avoid unnecessarily obscuring the present invention. 

[20] Referring now to the drawings, wherein like reference numerals designate identical or 
corresponding parts throughout the several views, and more particularly to Figure 1 thereof, 
there is illustrated a system 100 in which a demodulator with communications link adaptation 
according to the present invention may be employed. In Figure 1, in the system 100 
according to the present invention, a network operations control center 104 transmits 
information on satellite uplink channel 106, such as received from sources 102 (e.g., the 
Internet, an Intranet, content sources, etc.), to a satellite 108. The satellite 108 then transmits 
modulated information (e.g., using Quadrature Phase Shift Keying (QPSK), etc.) on a device 
downlink channel 1 10 to a device 112, such as a personal digital assistant (PDA), a personal 
information assistant (PIA), a personal computer (PC), a laptop PC, a television, an Internet 
appliance, a cellular phone, a set-top box, etc. 

[21] The device 112 includes an antenna 1 12a and a satellite communications transceiver 
(not shown) and thus is able to receive the modulated information on the downlink channel 
110. Such a satellite communications transceiver may include the demodulator with 
communications link adaptation according to the present invention, as will be further 
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described in detail with respect to Figure 2, to demodulate the information received on the 
downlink channel 110. 

[22 J The device 112 may make requests for information and/or transmit information via a 
device uplink channel 114. The satellite 118 receives information transmitted from the 
device 1 12 on the device uplink channel 1 14 and transmits the received information to the 
network operations control center 104 via a satellite downlink channel 116. The network 
operations control center 104 then may forward the information received on the satellite 
downlink channel 116 from the satellite 1 18 to the sources 102 (e.g., the Internet, an Intranet, 
content sources, etc.). 

[23] With the above-noted system 100, video download, audio download, graphics 
download, file download, pay per view, video-on-demand, audio-on-demand, Internet surfing, 
e-mail, voice communications, text communications,. paging functions, etc., maybe 
implemented on the device 112. One or more interface mechanisms may be used in the 
system 100, for example, including Internet access, telecommunications in any form (e.g., 
voice, modem, etc.), wireless communications media, etc., via the communication network 
104 and the satellite communications channels 106, 1 10, 1 14, and 1 16. The system 100 
information also may be transmitted via direct mail, hard copy, telephony, etc., when 
appropriate. 

[24 J Accordingly, the systems 104, 108 and 112 may include any suitable servers, 
workstations, personal computers (PCs), laptop PCs, personal digital assistants (PDAs), 
Internet appliances, set top boxes, other devices, etc., capable of performing the processes of 
the present invention. The systems 104, 108 and 112 may communicate with each other 
using any suitable protocol and, for example, via the communications network 102 and the 
communications channels 1 06, 1 10, 1 14 and 116 and may be implemented using the 
computer system 901 of Figure 9, for example. 

[25 J It is to be understood that the system in Figure 1 is for exemplary purposes only, as 
many variations of the specific hardware used to implement the present invention are 
possible, as will be appreciated by those skilled in the relevant art(s). For example, the 
functionality of the one or more of the systems 104 and 108 may be implemented via one or 
more programmed computers or devices. To implement such variations as well as other 
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variations, a single computer (e.g., the computer system 901 of Figure 9) may be programmed 
to perform the special purpose functions of, for example, the systems 104 and 108 shown in 
Figure 1 . On the other hand, two or more programmed computers or devices, for example as 
in shown Figure 9, may be substituted for any one of the systems 104, 108 and 1 12. 
Principles and advantages of distributed processing, such as redundancy, replication, etc., 
may also be implemented as desired to increase the robustness and performance of the system 
100, for example. 

126] The communications network 102 may be implemented via one or more 
communications networks (e.g., the Internet, an Intranet, a wireless communications network, 
a satellite communications network, a cellular communications network, a hybrid network, 
etc.), as will be appreciated by those skilled in the relevant art(s). In a preferred embodiment 
of the present invention, the communications network 102 and the communications channels 
106, 110, 114 and 116 and the systems 104, 108 and 112 preferably use electrical, 
electromagnetic, optical signals, etc., that carry digital data streams, as are further described 
with respect to Figure 9. The demodulator according to the present invention will now be 
described in detail in the following sections and with reference to Figures 1-9. 

[27J Figure 2 is a block diagram illustrating a demodulator 204 with communications link 
adaptation, which may be used in the system of Figure 1, according to the present invention. 
A matched filter 202 receives a signal over the communications channel 110 and outputs a 
filtered signal r{t). The signal r(t) filtered by the matched filter 202 is then passed to the 
demodulator 204 with communications link adaptation. The demodulator 204 with 
communications link adaptation receives the filtered signal r{f) and demodulates (e.g., using 
Bayesian demodulation techniques, etc.) the filtered signal r(f) to generate hard decisions 
used downstream, for example, for burst extraction, payload extraction, etc. The operation of 
the demodulator 204 of the present invention with communications link adaptation will now 
be described in detail. 

[28 J Generally, the demodulator 204 of the present invention may be employed in a 
wireless communications network, such as a satellite-based time-division-multiple-access 
(TDMA) GPRS wireless network of Figure 1 . It is assumed that there is only one path (i.e., 
Rician flat fading) in the satellite device downlink channel 110, which is quite common for 
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satellite communications. In such an application, for example, a Bayesian demodulation 
scheme, etc., may be used in the demodulator 204 in, for example, GPRS wireless 
communications. The demodulator 204 adaptively estimates the parameters for link 
adaptation. The demodulator 204 demodulates a received signal distorted by Rician fading 
based on, for example, one GPRS burst. An important parameter for link adaptation is signal 
to-noise ratio (SNR) known as signal quality indicator (SQI). Most modulator designs treat 
link adaptation and demodulation as two separate tasks. The present invention, however, 
unifies link adaptation and demodulation into one framework. A Gaussian distribution may 
approximate the physical channel of a satellite-based GPRS network. Then, according to the 
present invention, for example, a Bayesian classifier may be applied to the demodulator 204 
design based on such approximation. The signal and noise levels are adaptively"5stimated 
using unsupervised clustering, such as Fuzzy c-Means (FCM), etc., as will now be described. 

UNSUPERVISED CLUSTERING - FUZZY C-MEANS (FCM) 
[29] FCM clustering is a data clustering technique, wherein each data point belongs to a 
cluster to a degree specified by a membership grade. This technique was originally 
introduced by Bezdek [2] as an improvement on earlier clustering methods. FCM clustering 
will now be briefly summarized. 

Definition 1 (Fuzzy c-Partition): 

[30] Let X = x h x 2 , ...,x n be any finite set, V cn be the set of real c xn matrices, and c be an 
integer, where 2 <c<n. The Fuzzy c-partition space for X is the set: 



[31] The row i of matrix U <e M fc contains values of the ith membership function, u u in the 
fuzzy c-partition U of X. 



c 



n 



M fc = Us N cn \u ik € [0, 1] V/, k\ where 




(1) 
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Definition 2 (Fuzzy c-Means Functionals) [2]: 

132} Let J m : M /c *R cp -^R + be: 



n 



c 



(2) 



k=\ i=l 



where U e Mf C is a fuzzy c-partition o/X.; v = (vj, v^, . . \ c ) e R cp , where V/ e is f/*e 
cluster center of prototype uu 1 <i <c; 



where \[\\ is any inner product induced norm on R p ; weighting exponential m e [1, oo); Uik is 
the membership of x* in fuzzy cluster u t ; andJ m (U, v) represents the distance from any given 
data point to a cluster weighted by that point's membership grade. 

[33J The solutions of: 

min 

UeM„veR°> W ' V > (Z 

are least-squared error stationary points of J m . An infinite family of fuzzy clustering 
algorithms, one for each m e (1, oo), is obtained using the necessary conditions for solutions 
of equation (4), as summarized by the following Theorem: 



(^) 2 =|| x ,_ Vf || 2 



(3) 
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Theorem 1 [2]: 

[34] Assume ||*|| to be an inner product induced norm: fix m e (1, co), let X have at least c 
< n distinct points, and define the sets (Vk): 



h = {i\l<i<c;d ik =\\K k -Y i \\ = Q} (5) 
{l,2,...,c}-/* (6) 



Then (U, v) e xR cp is globally minimal for J m only if (wherein (f) denotes an empty set): 



or 



h^t <j>=$ uik= 0 V/ e 7* awd ^jT^ 1, 

iel k 



(8) 



2] (nor* / ^] («*r vi 



(9) 



*=1 
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[35] Bezdek proposed the following iterative method [2] to minimize J m (\J 9 v): 

1 . Fix c, 2 < c < n; choose any inner product nonn metric for R p ; and fix m, 1 < m < oo. 
Initialize l/ 0) e (e.g., choose its elements randomly from the values between 0 and 1). 
Then at step / (/ =1,2, ...): 

2. Calculate the c fuzzy cluster centers v,-® using equation (9) and U^. 

3. Update using equations (7) or (8). 

4. Compare to U (M) using a convenient matrix norm, i.e., if ||U (/) ~U (/_1) || < Si stop; 
otherwise, return to step 2. 

[36] The present invention may apply the FCM method to cluster, for example, one GPRS 
burst cell to 4 clusters. This is because, for example, Quadrature Phase Shift Keying (QPSK) 
modulation is used. The details are presented in the following section. However, the FCM 
method may be applied to other types of modulation techniques, such -as M-ary Phase-Shift 
Keying (M-PSK), Quadrature Amplitude Modulation (QAM), Pulse Amplitude Modulation 
(PAM), etc., as will be appreciated by those skilled in the relevant art(s). 
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SYSTEM MODEL 

[37] A satellite channel is often modeled as a Rician fading channel. Rician fading occurs 
when there is a strong specular (i.e., direct path or line of sight component) signal in addition 
to the scatter (i.e., multipath) components. The channel gain: 

g(t)=gi(t)+jg Q (t) (10) 

may be treated as a wide-sense stationary complex Gaussian random process, wherein g#) 
and gQ(t) may be Gaussian random processes with non-zero means mj{t) and rriQ(t), 
respectively. Because gj(t) and g Q (i) may have a same variance cr g 2 , the magnitude of the 
received complex envelope may have a Rician distribution [5] : 

x x + s 2 xs 

Pctx)=-— exp{~ — — }/o(— ) x>0 (11) 
o 2 2c? o 2 

where: 

s 2 = m 2 (t) + m Q 2 (t) (12) 
and /o(') may be a zero order modified Bessel function. 

[38] Such a channel may be known as a Rician fading channel. A Rician channel typically 
is characterized by two parameters, Rician factor K, which is the ratio of the direct path power 
to that of the multipath, i.e., K = s 2 /2c? [5] and the Doppier spread (or single-sided fading 
bandwidth) f d . The Rician fading may be simulated using a direct path added by a Rayleigh 
fading generator. The Rayleigh fade generator may be based on Jakes' model [3] in which an 
ensemble of sinusoidal waveforms are added together to simulate the coherent sum of 
scattered rays with Doppier spread^ arriving from different directions to the receiver. The 
amplitude of the Rayleigh fade generator is controlled by the Rician factor K. The number of 
oscillators to simulate the Rayleigh fading is, for example, 60. 
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[39] In the present invention, a system model 300, corresponding to the system 100 of 
Figure 1, as shown in Figure 3, may be used for simulation (e.g., using Mathcad by MathSoft 
Engineering & Education, Inc., etc.). In Figure 3, the system model for the transmitter (e.g., 
provided in the satellite 108) may include a random bits generator 302, a scalar-to-vector 
converter 304, a burst generator 306, a bits to integer converter 308, a modulator 310, an up- 
sampler (e.g., by 16) 312, a pulse shaping filter 314 (e.g., a square root raised cosine filter 
with roll off factor 0.35), and output of a Rician fading channel model 316 summed via 
summer 320 with output of a complex white noise generator 318. The system model for the 
receiver (e.g., provided in the device 112) may include the matched filter 202, a down- 
sampler 322 (e.g., by 16), a vector-to-scalar converter 324, the demodulator 204 (e.g., 
Bayesian, etc.), a burst extractor 326 and a bit error counter 328. 

[40} In Figure 4, a burst format used is shown. In Figure 4, the burst format may include, 
for example, 468 QPSK symbols per burst, 10 guard symbols at the beginning and end of the 
burst; 5 symbols including unique words (Uws, also referred to as a "training sequence," 
"synchronization words," etc.) for training; 24 public user information (PUT) symbols and 419 
symbols for payload. The random bits generator 302 generates, for example, a binary data 
stream with equally likely zeros and ones, which may be used for the payload bits (e.g. y 838 
bits). The burst builder 306 may insert some header and control bits and make a complete 
burst with, for example, 936 bits. The bits then are modulated to, for example, 468 QPSK 
symbols via the modulator 310. 
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LINK ADAPTATION AND BAYESIAN DEMODULATOR FOR GPRS 
Theoretical Basis 

[41 1 The matched filter 202 output when sampled in time-synchronization may be modeled 
as: 



r{i)=g(f)s{f) + n{t) (13) 



where: 



n(t) = m(t)+jndt) (14) 



is additive white Gaussian noise (AWGN) with mean of 0 and variance <j 2 in the in-phase (I) 
and quadrature (Q) components. For QPSK modulation, s(t) e {l 9 j 9 -1, -j} may be the 
signal points. Based on different values of s(t), the following results may be derived: 

[42] 1. KX0 = l,then: 



Hi) = g{i) + n(f) (15) 
= g&)+jg&t) + nAt)+jndi) (16) 
= [gM + nm+AgQ^ + riQit)] (17) 

[43] Since both gj{t) and n£t) may be Gaussian distributions with mean mj{t) and 0 and 
with variance a g 2 and a 2 , respectively, rj(t) A g#) + n^t) may be a Gaussian distribution with 
mean mXO and variance a g 2 + <j n 2 [1]. Similarly, r^(/) a + //g(£) may be a Gaussian 
distribution with mean mg(0 and variance c% 2 + a 2 . 
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[44] 2. If =y: 

= Jg(0 + n(t) (18) 
= + n0) +jn Q (t) (19) 

= [~gQ{t) + nit)]+jMt) + n^t)] (20) 

then rXO 4 -ggjit) + n0) may be a Gaussian distribution with mean -m^t) and variance cr/ + 
On and r^(?) A gXO + w e(0 ma Y De a Gaussian distribution with mean mX^) and variance o g 2 + 

O n \ 

145] 3. Ifs(0 = -1: 



K0 = -g(t) + n(t) (21) 
= -gtf)-jg(tf) + ntf)+JnQLf) ■ (22) 

= hg#) + »K0] + A-gQ(t) + »e(0] (23) 

then rXO 4 -g'XO + W X0 ma Y De a Gaussian distribution with mean -mi(t) and variance o g 2 + 

cr„ 2 and r e (f) 4 -gg(/) + «g(0 ma y De a Gaussian distribution with mean -mg(/)'and variance 

2 , 2 
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[46] 4. l£s(t) = -j: 

Ki) = -jg(t) + n(t) (24) 
= gQ{i) -jgtf) + nit) +jn Q {t) (25) 

= \ggff) + nit)} +j[-git) + jfgCO] (26) 

then rXO 4 g<2(/) + W X0 ma Y be a Gaussian distribution with mean m^t) and variance a g + 
ct^ 2 and r^(/) 4 -gXO + w f>(0 ma y ^ e a Gaussian distribution with mean -mif) and variance 

2 2 

[47] Summarizing the above results, for r(£) = rXO + j^Q{t\ equations (27) may be 
obtained, as follows: 

rit) ~ N{-mii), a 2 + o^ 2 ), r^) - N{-m^t\ cr 2 + o-„ 2 ), if - 1 

rXO ~ N{"~m Q {t\ cr 2 + a-, 2 ), r Q {t) ~ iV(-;mi(t) 5 cy g 2 +a n 2 ), if *(0 =7 " 

rXO - N(-;-mtf), cr g + <x„ 2 ), - N(-;-m Q (i), <r 2 + <t„ 2 ), if *(*) = -1 

rXt) ~ N(';m<#) 9 a 2 + <7„ 2 ), rg(0 - NQ'-m^t), a 2 + <x„ 2 ), if s(t) = -j 

where 7V(*; m, a 2 ) may denote a Gaussian distribution with mean m and variance a 2 . 
Accordingly, the received signals of one burst may be clustered to four clusters and the 
signals associated with each cluster may have a Gaussian distribution. The mean (i.e., time 
average) and variance of each cluster then may be used to compute the signal and noise level, 
i.e., to estimate the SQI and the demodulator 204 (e.g., Bayesian, etc.) maybe implemented 
based on this technique. 
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ESTIMATION OF THE SQI AND 
DESIGN OF THE BAYESIAN DEMODULATOR 
Determining the Mean, Variance, SQI, and Cluster Category 

[48] There may be, for example, 5 QPSK symbols that may be used to determine the 
category of each cluster. Suppose the 5 QPSK symbols are all 1 's. In this case, the FCM 
method may be used to cluster the 468 symbols (i.e., one burst) into c = 4 clusters, wherein 
each cluster has the mean (v h i = 1, 2, 3, 4). The FCM method may also generate U, a 4x468 
matrix in such application. Every received symbol r k (k= 1,2, 468) may have four 
membership grades u& e U (i = 1, 2, 3, 4 and ZVi Wjjk =1) corresponding to the four clusters. 
Based on the four values of u ik (i = 1, 2, 3, 4) for each k, which cluster each symbol belongs to 
maybe determined based on the maximum membership in u& (i = 1, 2, 3, 4). Accordingly, 
the 468 symbols may be clustered to 4 clusters using such a hard decision, which may be 
generated by the demodulator 204. 
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[49 j The variance of each cluster is computed based on such decision. The SQI may be 
estimated using: 



1 A 

Eb ~ E : 



S 



A 



2 H X (^ /2+CT -G 2 ) 



(28) 



A'" 4 

11 

V V (m in f+m in ^) 

2 4N , ri 

n=l i=l 

= lOlogjo ~ ~ (29) 

1 

-TV (o^+o-mg 2 ) 

v W=l 1=1 



TV 4 

^ ^ (tn in I 2 + m in (?) 
n=l i-1 

= lOlogio (30) 

TV 4 



where [m^Z, /w,-„Q], and [<7 in I, &inQ] may denote the estimated mean and standard deviation 
(std) of the I and Q components of the rth cluster in the nth burst, respectively. 
[50] Based on the cluster where the normalized UW signal has been assigned (i.e., based 
on the maximum membership), such a cluster may be assigned to the category "1." Because 
of the channel fading and randomness of the noise, normalized UW signals may be clustered 
to different clusters and majority logic may be used to determine which cluster may be 
assigned to the category "1." Once the cluster with category "1" is determined, the remaining 
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three clusters maybe assigned to the "j", and "-j" categories, for example, in a 
counterclockwise order from the cluster assigned category "1." 

Bayesian Demodulator Computation Formula 

[51] Bayesian detection may then be applied to every signal point in the burst using the 
following rule: 

[52] The signal r(t) may given by a t (i = 1, 2, 3, 4), where a t e { 1 ,y, - 1 , -j) if: 



p(r(t)\s(t) = at) >p(r(t)\s(t) = aj) Vaj * a t (31) 
[53J To compute p{r{i)\s{i) = at), r 4 [rj(t), r Q (t)Y may be given by: 



p{r{t)\s(t) = ai) =p(r\ad 
(2tc)|E ; | 1/2 eXpt ~ 2 (r ~ mi ) TE ' _1 ( r 1 m ')J 



(32) 
(33) 



where m, 4 [m/, mfiy and E/ = diag{a g 2 + a n 2 , cr g 2 + a 2 } are the mean vector (2x 1) and 
covariance matrix (2x2) of [rfc), r e (r)] r . 

[54] Figure 5 is a flow chart illustrating the operation of the demodulator 204 with 
communications link adaptation of Figure 2, according to the present invention. In Figure 5, 
at step 502, the demodulator 204 receives the signal r(t) from the matched filter 202 based on 
the signal received over the device downlink channel 1 10 by the matched filter 202. At step 
504, clusters are extracted from the received signal r(t) using unsupervised clustering, such as 
FCM, etc., as previously described. At step 506, the mean and standard deviation are 
computed for each extracted cluster, as previously described. At step 508, a category for each 
extracted cluster (e.g., "1," "j", "-1", and "-j" for QPSK modulation) is determined, as 
previously described. At step 510, demodulation (e.g., Bayesian, etc.) is performed based on 
the computed mean, standard deviation and determined category to generate the hard 
decisions, as previously described, completing the operation. 
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SIMULATIONS 

[55] For simulation purposes, a Rician fading channel with, for example, a Rician factor K 
= 9dB, a Doppler shift f d = 10Hz and symbol rate 93.6ks/s (i.e., 5ms per burst) may be used. 
The information (i.e., payload) bit rate of \61.6kb/s may be used. How well the FCM 
algorithm works in estimating the mean and variance for the four clusters (i.e., in QPSK 
modulation) may be evaluated in terms of root-mean-square-error (RMSE), which may be 
defined as: 




N 4 
] A 



RMSE m = *\ I — Z Z K w ^-»^ + ( m ™Q -™;nQf] (34) 
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N 4 



RMSE std = *\ / — Z Z ( ^ / -^ 7 > 2 + ( <**G -°inQfl (35) 



47^ 



where [/w'V, m A y/2 <2] and [a A z >/, o A in Q] denote the estimated mean and standard deviation 
(std) of the I and Q components of the zth cluster in the nth cell respectively; [m in I, m in Q] and 
[<j in I, <SinQ\ are the actual mean and std of the I and Q components (i.e., obtained via 
supervised clustering in which we typically know the exact category of received signal) of the 
ith cluster in the nth cell, respectively; and N is the total number of cells. 
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[561 Figures 6a and 6b are graphs showing the performance (i.e., RMSE m and RMSE std ) of 
FCM clustering for 5,000 cells at each EM (dB) value, which may be defined as: 
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[57] The simulation system 300 may be calibrated and the average signal power E s may be 
normalized to 1 before the demodulator 204. From Figures 6a and 6b, it is observed that both 
RMSE m and RMSE std are quite small, which means that the mean and variance of the received 
symbols obtained via the FCM approaches those obtained via supervised clustering. 
Accordingly, when the symbol categories of received signals are unknown, FCM clustering 
may be used to extract the mean and variance of four clusters, which makes the optimal 
(Bayesian) demodulation possible. 

| 58 j In Figure 7, a plot of the average SQI estimation error (in dB) versus the actual Ei/No 
(dB) for 5,000 bursts is shown. From Figure 7, it is observed that the error is quite small 
especially for high Ef/No, which means that our FCM-based clustering can be used to estimate 
the SQI and evaluate the link quality. 

[59] The FCM-based demodulator 204 of the present invention using, for example, 
Bayesian demodulation was compared against a block phase estimation (BPE)-based 
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demodulator [6] and the results are shown in Figure 8. BPE is a semi-blind demodulator, 
which uses a maximum-likelihood method to estimate a phase gain of a channel and then uses 
a known sequence (e.g., 5 QPSK unique words) to remove phase ambiguity because of an 
inverse-tangent function used in the phase estimation. This method has been widely used in 
mobile communications with burst transmission. 

[60] For such a channel, simulations for different EM values were performed. At each 
EiJNq value, the simulations for 5000 bursts were ran and the average bit error rate (BER) for 
the FCM-based Bayesian demodulator and BPE-based demodulator were obtained. The 
performances are plotted in Figure 8. Also plotted was the theoretical BER at K= 9dB. From 
Figure 8, it is observed that the FCM-based demodulator 204 of the present invention 
performs better than the BPE-based demodulator, achieving a gain of 0.2dB andlipproaching 
the theoretical BER. 

[61 ] The present invention stores information relating to various processes described 
herein. This information is stored in one or more memories, such as a hard disk, optical disk, 
magneto-optical disk, RAM, etc. One or more databases, such as the databases within the 
systems 104, 108 and 1 12, etc., may store the information used to implement the present 
invention. The databases are organized using data structures (e.g., records, tables, arrays, 
fields, graphs, trees, and/or lists) contained in one or more memories, such as the memories 
listed above or any of the storage devices listed below in the discussion of Figure 9, for 
example. 

[62] The previously described processes include appropriate data structures for storing data 
collected and/or generated by the processes of the system 100 of Figure 1 in one or more 
databases thereof. Such data structures accordingly will includes fields for storing such 
collected and/or generated data. In a database management system, data is stored in one or 
more data containers, each container contains records, and the data within each record is 
organized into one or more fields. In relational database systems, the data containers are 
referred to as tables, the records are referred to as rows, and the fields are referred to as 
columns. In object-oriented databases, the data containers are referred to as object classes, 
the records are referred to as objects, and the fields are referred to as attributes. Other 
database architectures may use other terminology. Systems that implement the present 
invention are not limited to any particular type of data container or database architecture. 
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However, for the purpose of explanation, the terminology and examples used herein shall be 
that typically associated with relational databases. Thus, the terms "table," "row," and 
"column" shall be used herein to refer respectively to the data container, record, and field. 

[63] The present invention (e.g., as described with respect to Figures 1-8) may be 
implemented by the preparation of application-specific integrated circuits or by 
interconnecting an appropriate network of conventional component circuits, as will be 
appreciated by those skilled in the electrical art(s). In addition, all or a portion of the 
invention (e.g., as described with respect to Figures 1-8) maybe conveniently implemented 
using one or more conventional general purpose computers, microprocessors, digital signal 
processors, micro-controllers, etc., programmed according to the teachings of the present 
invention (e.g., using the computer system of Figure 9), as will be appreciated by those skilled 
in the computer and software art(s). Appropriate software can be readily prepared by 
programmers of ordinary skill based on the teachings of the present disclosure, as will be 
appreciated by those skilled in the software art. Further, the present invention may be 
implemented on the World Wide Web (e.g., using the computer system of Figure 9). 

[64] Figure 9 illustrates a computer system 901 upon which the present invention (e.g., 
systems 104, 108, 112, etc.) can be implemented. The present invention may be implemented 
on a single such computer system, or a collection of multiple such computer systems. The 
computer system 901 includes a bus 902 or other communication mechanism for 
communicating information, and a processor 903 coupled to the bus 902 for processing the 
information. The computer system 901 also includes a main memory 904, such as a random 
access memory (RAM), other dynamic storage device (e.g., dynamic RAM (DRAM), static 
RAM (SRAM), synchronous DRAM (SDRAM)), etc., coupled to the bus 902 for storing 
information and instructions to be executed by the processor 903. In addition, the main 
memory 904 can also be used for storing temporary variables or other intermediate 
information during the execution of instructions by the processor 903. The computer system 
901 further includes a read only memory (ROM) 905 or other static storage device (e.g., 
programmable ROM (PROM), erasable PROM (EPROM), electrically erasable PROM 
(EEPROM), etc.) coupled to the bus 902 for storing static information and instructions. 
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[65] The computer system 901 also includes a disk controller 906 coupled to the bus 902 to 
control one or more storage devices for storing information and instructions, such as a 
magnetic hard disk 907, and a removable media drive 908 (e.g., floppy disk drive, read-only 
compact disc drive, read/write compact disc drive, compact disc jukebox, tape drive, and 
removable magneto-optical drive). The storage devices may be added to the computer system 

901 using an appropriate device interface (e.g., small computer system interface (SCSI), 
integrated device electronics (IDE), enhanced-IDE (E-IDE), direct memory access (DMA), or 
ultra-DMA). 

[66] The computer system 901 may also include special purpose logic devices 918, such as 
application specific integrated circuits (ASICs), full custom chips, configurable logic devices 
(e.g., simple programmable logic devices (SPLDs), complex programmable logic devices 
(CPLDs), field programmable gate arrays (FPGAs), etc.), etc., for performing special 
processing functions, such as signal processing, image processing, speech processing, voice 
recognition, infrared (ER) data communications, satellite communications transceiver 
functions, the demodulator 204 functions, etc. 

[67] The computer system 901 may also include a display controller 909 coupled to the bus 

902 to control a display 910, such as a cathode ray tube (CRT), liquid crystal display (LCD), 
active matrix display, plasma display, touch display, etc., for displaying or conveying 
information to a computer user. The computer system includes input devices, such as a 
keyboard 911 including alphanumeric and other keys and a pointing device 912, for 
interacting with a computer user and providing information to the processor 903. The 
pointing device 912, for example, may be a mouse, a trackball, a pointing stick, etc., or voice 
recognition processor, etc., for communicating direction information and command selections, 
to the processor 903 and for controlling cursor movement on the display 910. In addition, a 
printer may provide printed listings of the data structures/information of the system shown in 
Figures 1-8, or any other data stored and/or generated by the computer system 901 . 

[68] The computer system 901 performs a portion or all of the processing steps of the 
invention in response to the processor 903 executing one or more sequences of one or more 
instructions contained in a memory, such as the main memory 904. Such instructions may be 
read into the main memory 904 from another computer readable medium, such as a hard disk 
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907 or a removable media drive 908. Execution of the arrangement of instructions contained 
in the main memory 904 causes the processor 903 to perform the process steps described 
herein. One or more processors in a multi-processing arrangement may also be employed to 
execute the sequences of instructions contained in main memory 904. In alternative 
embodiments, hard- wired circuitry may be used in place of or in combination with software 
instructions. Thus, embodiments are not limited to any specific combination of hardware 
circuitry and software. 

[69] Stored on any one or on a combination of computer readable media, the present 
invention includes software for controlling the computer system 901, for driving a device or 
devices for implementing the invention, and for enabling the computer system 901 to interact 
with a human user (e.g., a user of the systems 104, 108, 1 12, etc.). Such software may 
include, but is not limited to, device drivers, operating systems, development tools, and 
applications software. Such computer readable media further includes the computer program 
product of the present invention for performing all or a portion (if processing is distributed) of 
the processing performed in implementing the invention. Computer code devices of the 
present invention may be any interpretable or executable code mechanism, including but not 
limited to scripts, interpretable programs, dynamic link libraries (DLLs), Java classes and 
applets, complete executable programs, Common Object Request Broker Architecture 
(CORBA) objects, etc. Moreover, parts of the processing of the present invention maybe 
distributed for better performance, reliability, and/or cost. 

[70] The computer system 901 also includes a communication interface 913 coupled to the 
bus 902. The communication interface 913 provides a two-way data communication coupling 
to a network link 914 that is connected to, for example, a local area network (LAN) 915, or to . 
another communications network 916 such as the Internet. For example, the communication 
interface 913 may be a digital subscriber line (DSL) card or modem, an integrated services 
digital network (ISDN) card, a cable modem, a telephone modem, etc., to provide a data 
communication connection to a corresponding type of telephone line. As another example, 
communication interface 913 may be a local area network (LAN) card (e.g., for Ethernet™, 
an Asynchronous Transfer Model (ATM) network, etc.), etc., to provide a data 
communication connection to a compatible LAN. Wireless links can also be implemented. 
In any such implementation, communication interface 913 sends and receives electrical, 
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electromagnetic, or optical signals that carry digital data streams representing various types of 
information. Further, the communication interface 913 can include peripheral interface 
devices, such as a Universal Serial Bus (USB) interface, a PCMCIA (Personal Computer 
Memory Card International Association) interface, etc. 

[71 1 The network link 914 typically provides data communication through one or more 
networks to other data devices. For example, the network link 914 may provide a connection 
through local area network (LAN) 915 to a host computer 917, which has connectivity to a 
network 916 (e.g. a wide area network (WAN) or the global packet data communication 
network now commonly referred to as the "Internet") or to data equipment operated by 
service provider. The local network 915 and network 916 both use electrical, 
electromagnetic, or optical signals to convey information and instructions. The signals 
through the various networks and the signals on network link 914 and through 
communication interface 913, which communicate digital data with computer system 901, are 
exemplary forms of carrier waves bearing the information and instructions. 

[72] The computer system 901 can send messages and receive data, including program 
code, through the network(s), network link 914, and communication interface 913. In the 
Internet example, a server (not shown) might transmit requested code belonging an 
application program for implementing an embodiment of the present invention through the 
network 916, LAN 915 and communication interface 913. The processor 903 may execute 
the transmitted code while being received and/or store the code in storage devices 907 or 908, 
or other non-volatile storage for later execution. In this manner, computer system 901 may 
obtain application code in the form of a carrier wave. With the system of Figure 9, the 
present invention may be implemented on the Internet as a Web Server 901 performing one or. 
more of the processes according to the present invention for one or more computers coupled 
to the Web server 901 through the network 916 coupled to the network link 914. 

[73] The term "computer readable medium" as used herein refers to any medium that 
participates in providing instructions to the processor 903 for execution. Such a medium may 
take many forms, including but not limited to, non- volatile media, volatile media, 
transmission media, etc. Non-volatile media include, for example, optical or magnetic disks, 
magneto-optical disks, etc., such as the hard disk 907 or the removable media drive 908. 
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Volatile media include dynamic memory, etc., such as the main memory 904. Transmission 
media include coaxial cables, copper wire, fiber optics, including the wires that make up the 
bus 902. Transmission media can also take the form of acoustic, optical, or electromagnetic 
waves, such as those generated during radio frequency (RF) and infrared (IR) data 
communications. As stated above, the computer system 901 includes at least one computer 
readable medium or memory for holding instructions programmed according to the teachings 
of the invention and for containing data structures, tables, records, or other data described 
herein. Common forms of computer-readable media include, for example, a floppy disk, a 
flexible disk, hard disk, magnetic tape, any other magnetic medium, a CD-ROM, CDRW, 
DVD, any other optical medium, punch cards, paper tape, optical mark sheets, any other 
physical medium with patterns of holes or other optically recognizable indicia, aRAM, a 
PROM, and EPROM, a FLASH-EPROM, any other memory chip or cartridge, a carrier wave, 
or any other medium from which a computer can read. 

[74] Various forms of computer-readable media may be involved in providing instructions 
to a processor for execution. For example, the instructions for carrying out at least part of the 
present invention may initially be borne on a magnetic disk of a'remote computer connected 
to either of networks 915 and 916. In such a scenario, the remote computer loads the 
instructions into main memory and sends the instructions, for example, over a telephone line 
using a modem. A modem of a local computer system receives the data on the telephone line 
and uses an infrared transmitter to convert the data to an infrared signal and transmit the 
infrared signal to a portable computing device, such as a personal digital assistant (PDA), a 
laptop, an Internet appliance, etc. An infrared detector on the portable computing device 
receives the information and instructions borne by the infrared signal and places the data on a 
bus. The bus conveys the data to main memory, from which a processor retrieves and 
executes the instructions. The instructions received by main memory may optionally be 
stored on storage device either before or after execution by processor. 

[75] Recapitulating, the present invention, advantageously, provides a demodulator 204 
(e.g., using Bayesian demodulation, etc.) including a link adaptation scheme that maybe 
employed, for example, in satellite-based General Packet Radio Services (GPRS) wireless 
network. The received signal of a satellite-based GPRS channel may be approximated by 
Gaussian distributions and then the demodulator 204 based on such distributions may be 
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implemented. The average signal and noise powers in one GPRS burst may be estimated 
using an unsupervised clustering method, such as fuzzy c-means (FCM), etc. Based on the 
estimated signal level and noise standard deviation, a signal quality indicator (SQI) may be 
evaluated and the parameters associated with the demodulator 204 may be implemented. The 
demodulation scheme uses link adaptation results, which simplifies the system design 
tremendously, as compared to conventional demodulators. Simulation results show that the 
link adaptation scheme works well and that the demodulator 204 of the present invention 
performs better than a block phase estimation (BPE) demodulator, achieving a gain of O.ldB. 
[76| The demodulator 204 of the present invention may be employed, for example, in an 
Inmarsat [9] multimedia communications system. Since the Bayesian demodulator typically 
| obtains a 0.2dB gain over the existing BPE demodulator, it can save millions of dollars in the 

satellite communications costs. The demodulator described in [8] is a training-based 
\ demodulator, wherein the number of unique words typically cannot be too small and the 
unique words typically may be located at several different locations (i.e., not in one place) in 
a burst so that the channel characteristics may be captured. The demodulator 204 of the 
present invention typically does not suffer from such constraints. 

[ 771 Although the present invention is described in terms of a demodulator used in a 
Quadrature Phase Shift Keying (QPSK) modulation environment, the present invention is 
applicable to other modulation environments, such as M-ary Phase-Shift Keying (M-PSK), 
Quadrature Amplitude Modulation (QAM), Pulse Amplitude Modulation (PAM), etc., as will 
be appreciated by those skilled in the relevant art(s). 

[ 78 J Although the present invention is described in terms of a demodulator used in a 
system using a satellite communications channel, the present invention is applicable to other 
systems that may employ a demodulator using other communications channels, such as a 
digital video broadcasting (DVB) communications channel, a terrestrial broadcast 
communications channel, a cellular communications channel, a Quadrature Phase Shift 
Keying (QPSK) communications channel, an M-ary Phase-Shift Keying (M-PSK) 
communications channel, a Quadrature Amplitude Modulation (QAM) communications 
channel, a Pulse Amplitude Modulation (PAM) communications channel, etc., as will be 
appreciated by those skilled in the relevant art(s). 
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[79] Although the present invention is described in terms of a demodulator using Bayesian 
demodulation techniques, the present invention is applicable to other types of demodulation 
techniques, as will be appreciated by those skilled in the relevant art(s). 

[80] While the present invention has been described in connection with a number of 
embodiments and implementations, the present invention is not so limited but rather covers 
various modifications and equivalent arrangements, which fall within the purview of the 
appended claims. 
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